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Superconducting Nb thin films with rectangular arrays of submicron antidots have been system- 
ically investigated by transport measurements. In low fields, the magnetoresistance curves demon- 
strate well-defined dips at integral and rational numbers of flux quanta per unit cell, which corre- 
sponds to a superconducting wire network-like regime. When the magnetic field is higher than a 
saturation field, interstitial vortices interrupt the collective oscillation in low fields and form vortex 
sublattice, where a larger magnetic field interval is observed. In higher fields, a crossover behavior 
from the interstitial sublattice state to a single-loop-like state is observed, characterized by oscilla- 
tions with a period of $o/""'"e//, originating from the existence of edge superconducting states with 
a size r e // around the antidots. 
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PACS numbers: 74.25. F-, 74.78.Na, 74.81.Fa 



I. INTRODUCTION 

Experiments on mesoscopic superconductors with di- 
mensions comparable to the superconducting charac- 
teristic length scales have demonstrated that the sam- 
ple topology strongly influences the superconducting 
properties, such as the phase boundary T C (H), magne- 
toresistance R(H), and field dependent critical current 
I C (H)^ ~— Various topologies (single loops,^ multiloops,— 
large infinite networks,— and arrays of antidots^) have 
been studied both experimentally and theoretically. 

In a perpendicular magnetic field, localized supercon- 
ducting state can first nucleate near the edge of samples 
within a thin layer of width W s ~ £,(T), similar to the nu- 
cleation of the surface superconductivity^ The so-called 
edge superconducting stat o 9 ' 10 has an enhanced critical 
field H C i{T). The enhancement of H C ^(T) above the bulk 
critical field H C 2(T) greatly depends on the curvature of 
the superconducting/normal interface and the surface-to- 
volume ratio £ii Much higher enhancement of the ratio 
H&l H C 2 up to 3.6 has been observed in Pb thin film with 
a dense square antidot latticed 

For the square arrays of antidots, when the narrowest 
separation AW between neighboring holes is smaller than 
a critical value 1.84£(T)j& nucleation is dominated by the 
thin wire-like edge supercondcting states and the cou- 
pling between them. This kind of arrays is well described 
by the theory of supercondcting wire networks! 12 ' 13 The 
oscillations of T C {H) or R(H) in such array are known as 
collective oscillations or network-like oscillations, whose 
period is corresponding to the area of the unit cell£ As 
shown by Bezryadin and Pannetier^ when the magnetic 
field is high enough, a crossover behavior in the T C {H) 
curve from collective oscillations to single-loop-like ('sin- 
gle object') oscillations has been observed. The period of 
the singlc-loop-likc oscillations is determined by the ef- 
fective hole (antidot) radius r e ff = rh + hW s , where rh is 



the hole radius. Due to the decrease of W s with increas- 
ing magnetic field, the single-loop-like oscillations are 
non-periodic. Although the dinrensional crossover behav- 
iors have been observed in square array of antidotsr^ 9 - 
few experiments on other symmetry of arrays are per- 
formed to check the universality. It is therefore inter- 
esting to study the rectangular arrays by introducing an 
anisotropy to the square arrays. We increase the length of 
one side of the rectangular unit cell, while keep the other 
one constant. In this way, we obtain a series of rectangu- 
lar arrays of antidots characterized by two features: (1) 
separation AW of the short side a is smaller than £(T); 
(2) AW of the long side b is larger than 4f (T) (i.e. little 
overlap of the edge states in this direction). Furthermore, 
for perforated samples, Abrikosov (interstitial) vortices 
can appear in the wide superconducting strips between 
holes, which are observed by imaging techniques^ ' 10 ' 14 ' 15 
Thus, for the rectangular arrays of antidots, an extra in- 
terstitial vortex state will be involved as magnetic field 
increases, in contrast to the crossover behavior found in 
the square arrays of antidotsi&i' 9 - 

In this paper, we have performed detailed systematic 
transport measurements on the rectangular arrays of an- 
tidots with various geometric parameters. We observe 
clearly successive crossover behaviors from network-like 
state to interstitial vortex state, then to the single-loop- 
like state. Hysteresis effect is observed for the interstitial 
vortex state. The crossover fields arc found strongly de- 
pendent on temperature, the aspect ratio of the unit cell, 
and hole size. 



II. EXPERIMENT 

High quality Nb thin films with a thickness of about 
100 nm were deposited on Si02 (300 nm)/Si substrates 
by magnetron sputtering. The Nb thin films have a crit- 
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ical temperature T c of 8.910 K and a superconducting 
transition width of about 21 mK (10% - 90%i? n crite- 
rion, where R n is the normal state resistance at 9 K). 
For standard transport measurements, three four-probe 
microbridges were patterned on a single Nb substrate 
with ultraviolet photolithography and etched by reac- 
tive ion etching (RIE) in O2 and SFg plasmas. Each 
Nb microbridgc had a 60 (im width for current flow- 
ing and a 60 /jm distance between two voltage connects. 
The substrates were spun with polymethyl metacrylate 
(PMMA) resist layer and baked at 170 °C for 1 minute. 
Designed patterns were written by electron beam lithog- 
raphy (EBL) on PMMA resist and developed for 40 sec- 
onds in MIBK:IPA (1:3) solution. Finally, the samples 
were etched by RIE and unexposed PMMA was removed 
in acetone. Rectangular arrays of circular antidots were 
obtained in the centers of the Nb microbridges. Note that 
one of the Nb microbridges was intentionally unexposed 
and used as a reference. 

Figure [TJ shows images of scanning electron microscope 
(SEM) and atomic force microscope (AFM) for sample 
SI, with a rectangular unit cell of 800 nm x 1200 nm. 
The narrowest separation AW between neighboring holes 
along the X direction is 50 nm (±4 nm). The images 
demonstrate that the overall periodicity is maintained 
very well. Sharp edges arc obtained after etching. 
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Figure 1. (color online) SEM and AFM images of a rectangu- 
lar antidot lattice(sample SI), with a unit cell of a x 6 (800 
nm x 1200 nm). The hole(antidot) radius is 375 nm and the 
distance between centers of holes is 800 nm, resulting in a 
width AW of the narrowest part of the constriction of 50 nm. 
The dashed rectangle indicates a unit cell. 

The measurements have been performed in Physical 
Properties Measurement System (PPMS-14, Quantum 
Design Inc.). Two phase lock-in amplifiers (SR830) are 
used for ac currents applied at a frequency of 30.9 Hz. 
The current is parallel to the long side of rectangular unit 
cell (Y direction). Applied magnetic field is perpendicu- 
lar to the film surface. We sweep the magnetic field with 
a step of 0.4 Oe in low field regime and 1 Oe in high field 
regime. The temperature stability is better than 2 mK 
during the measurements. The superconducting coher- 
ence length £(0) is 11.3 nm and the penetration depth 
A(0) is 74.0 nm, determined by measuring the T C (H) of 
the reference Nb microbridge.— Thus, we have £(£) = 113 
nm, A e //(t) ra \ 2 (t)/d = 1389.7 nm, at t = 0.990, where 
t = T/T c is the reduced temperature, T = 8.821 K, and 
d = 100 nm is the thickness of thin film. 
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Figure 2. (color online) (a) Magnetoresistance of sample SI 
is measured at T = 8.640 K (t = 0.993) with I = 30 fj,A. The 
dashed vertical lines divide R(H) curve into three regions. 
The downward arrows indicate the positions of the crossover 
fields Hs and Hd- Inset of (a): Enlarged plot of the low field 
regime. Fractional reduced magnetic fluxes are indicated with 
upward arrows, (b) Index of the R(T) peaks N as a function 
of the field position of the minima. Solid lines are linear fits 
for the data, with slopes (fcl, fc2 and fc3), corresponding to 
the three regions in upper panel. 



III. RESULTS AND DISCUSSION 
A. Magnetoresistance and hysteresis effect 

The field dependent R{H) curve for sample SI (800 x 
1200, r h = 375 nm) at T = 8.640 K and current I = 30 
^A is given in Fig. 12(a). The T c of sample SI is 8.702 K 
(50%i? n criterion). Three different regimes of the R{H) 
curve can be distinguished by the shape of dips(minima) 
and the magnetic field interval AH between two consec- 
utive dips: the low field (region-I), the intermediate field 
(region-II) and the high field (region-Ill) regimes. The 
downward arrows indicate the crossover fields at Hs = 
125 Oc and Ho = 352 Oe. The number N is a sequence 
number of dips relative to the one at zero field, which 
marks as N = 0. We plot the number N as a func- 
tion of the position of the dips, as shown in Fig. [2(b). 
The data are fitted by straight lines with different slopes 
(fcl = 1/20.8, fc2 = 1/34.9 and fc3 = 1/37.9), correspond- 
ing to the three regions in the upper panel. The slope 
fc is a reciprocal of the average field interval AH. It is 
found that AH in the three regions increases with field. 

In the low field regime (H < H s ), AH is 20.8 Oc, 



corresponding to one flux quantum per unit cell. It 
slightly deviates from the theoretic value Hr = $ /ab = 
21.6 Oe, where $0 = 20.7 G-/im 2 is the flux quantum.— 
The samples behave like the rectangular superconduct- 
ing wire network, owing to the separations of neighbor- 
ing holes comparable with £(T) near T c &H-~— The in- 
set of Fig. [Ha) shows a magnification of the low field 
R(H) curve. The fractional reduced magnetic flux / = 
$/$o = (1/4, 1/2, 3/4) is clearly visible. These fea- 
tures reflect the collective behavior of multiconnected su- 
perconducting wire network, which can be described by 
the mean field Ginzburg-Landau theory*^ 9 - The descrip- 
tion of the collective behavior is different from the arrays 
with smaller and weaker pinning centers. In the latter 
case, the observed oscillatory dips in the magnetoresis- 
tance were explained by the vortex (multiquanta- vortex) 
matching model in the London limit.— ~— 

The maxima of magnctorcsistance show approximately 
the same magnitude in the low field regime, suggesting 
that the multiquanta vortex is effectively confined in each 
large antidote When the magnetic field is larger than the 
saturation field Hs, the R(H) behavior changes drasti- 
cally. The collective oscillations are interrupted by the 
formation of additional vortices in the interstitial regions. 
The interstitial vortices cause phase decoherence events 
(i.e. phase slips^) or an incoherent 'glassy' configuration 
of order parameter phased resulting in broad dips, miss- 
ing of fine fractional structures, and increasing of field 
intervals. Similar phenomena have also been observed in 
rectangular arrays of magnetic dots, but explained by the 
reconfiguration transitions of vortex lattice! 21 ' 24 

In the high field regime (H > Ho ) , a rapid increase of 
resistance background is observed and AH = 37.9 Oe is 
larger than that in region-II. As we will discuss below, the 
large period oscillations in region-Ill correspond to the 
single-loop-like behavior. Thus, for the single-loop-like 
oscillations, a radius r$ deduced from (^q/ttAH) 1 ^ 2 = 
AY! nm is comparable to the one given by the expression 
of r e ff = (r h + AW) = 375 + 50 = 425 nm. It should be 
mentioned that we have replaced W s by AW, due to the 
fact that AW < £(T). Since 2r > a, it demonstrates 
that the supercurrents strongly overlap in the thin wires 
along the Y direction! 28 ' 29 To further study the R(H) 
behavior, we performed hysteresis measurements at the 
same temperature and current. It is expected that a 
hysteresis effect occurs in the intermediate regime owing 
to the existence of interstitial vortices. 

In Fig. [31 hysteresis curves are recorded in the follow- 
ing way: the magnetic field first increases from -100 Oe 
to 750 Oe, then decreases from 750 Oe to -100 Oc. The 
curves are reversible in region I and III. In rcgion-I, since 
only coreless vortices are present inside the holes, increas- 
ing or decreasing magnetic filed would produce the same 
magnetoresistance oscillations. In region-II, at first, ev- 
ery hole encloses 6 < £>o- There is a competition between in- 
creasing the flux per hole and accommodating vortices at 
interstitial positions in the wide strips. The whole system 
favors to stay in the lowest free energy state. Thus, the 
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Figure 3. (color online) Hysteretic effect of the magnetore- 
sistance for sample SI. Curves with open circles and filled 
squares correspond to increasing and decreasing field, respec- 
tively. The downward arrows mark the positions of Hs and 
Hd, which divide the curves into three parts: (I) H < Hs, 
network-like region, where no vortices are located outside the 
holes; (II) Hs < H < Hd, interstitial vortex sublattice ap- 
pears in the wide strips, with a normal core size of 2£(T); 
(III) H > Hd, the edge superconducting states are localized 
around each hole. The insets show schematic drawings of the 
vortex patterns for each region. 



additional vortices arc expected to appear in the weak 
superconducting regions. Furthermore, the interstitial 
vortices in the wide strips can form stable sublatticcs, 
which produce dips in R(H) curves . 22 ' 23 The increase 
in AH up to 34.9 Oc illustrates the presence of a com- 
pressed vortex distribution, compared to region-I. This 
increase is a little higher than the theoretical expecta- 
tion for reconfiguration of square lattice AH$ q = $o/a 2 
— 32.3 ± 1.7 Oei 21 ' 24 The hysteresis effect in region- 
II is due to the intrinsic pinning forces and barrier for 
vortex motion coming from remaining superconducting 
regions £2. With increasing field, the vortex patterns be- 
come more complex, and parts of the interstitial regions 
turn into normal state due to the penetration of the mag- 
netic field. In region-Ill, most of the regions in SI turns 
to normal states, except that edge states around holes 
are still superconducting. The curves become reversible 
again. A synchronized entrance of an additional vortex 
in each hole causes a resistance minimum in the R(H), 
like the Little-Parks effect* 1 - 

From the hysteresis measurements, the boundaries of 
three regions can be clearly distinguished. Similar results 
have also been observed on sample S2 (a = 800 nm, b = 
2000 nm, r h = 373 nm), with H s = 39 Oc and H D = 319 
Oe at T = 8.640 K. In contrast to the work on magnetic 
dots^ the hysteresis curves were measured just above 
or under the transition field, then the divisions of three 
regimes were possibly missed in their discussion. 



B. Parameters affecting the transitions 

Figure @] shows the R(f) curves of samples SI and S2 
measured at several temperatures, with / = 30 fiA [panel 
(a)] and 200 //A [panel (b)], respectively. Since these two 
samples show very similar temperature dependence of the 
R(f) curves, we focus our discussion on sample SI. In 
Fig.QJa), from top to bottom, the temperature decreases 
from 8.723 K to 8.540 K (partly shown in the figures). In 
the network-like region of R(f) curves, integer dips are 
always visible at these temperatures. Even some dips at 
fractional / = (1/4, 1/3, 1/2) fields are well developed 
when the temperature is lower than 8.680 K. However, 
the magnetoresistance oscillations in the intermediate re- 
gion and single-loop-like region are very sensitive to the 
temperature variation. At high temperatures close to T c 
(T > 8.680 K), the oscillations in these two regions are 
broad and shallow. With decreasing temperature, the 
oscillations become more pronounced. With further de- 
crease of the temperature, the oscillations become weaker 
and finally disappear. This illustrates different oscillation 
nature in these three regions. 

The saturation number Ng is defined by the largest 
possible number of the flux quanta trapped by an antidot. 
Mkrtchyan and Shmidlj 3 -^ have theoretically estimated 
the maximum possible number of vortices trapped by a 
single insulating inclusion with an expression of Ng t (t) = 
rfc/2£(£). In our case, larger hole and smaller separation 
of antidots along the X direction are both used, where 
the edge superconductivity is more significant and pro- 
duces effectively constraint to the flux^^ As indicated 
by the arrow in Fig. H^a), N s = 6 at T = 8.615 K. It 
is larger than the theoretical value Ng t (t = 0.99, r^ = 
375 Dm) ~ 2. Ng as a function of temperature for these 
two samples is shown in Fig. Eta), compared with the 
theoretical value Nst{t). As the temperature decreases, 
Ng increases stepwisely. Interestingly, the experimental 
result for sample S2 is more close to the theoretical value 
at low temperature (Ng rs 3, at t = 0.978). This is due 
to the fact that the distributions of antidots in sample 
S2 are more sparse along the Y direction (or a larger as- 
pect ratio b/a). In contrast, Ng for sample SI is nearly 
2^3 times of that for S2 at the same reduced tempera- 
ture. In an array of antidots, the antidot-vortex interac- 
tion strongly affects the saturation number, resulting in a 
larger saturation number in the dense arra y 32 ' 33 Besides 
that, the hole size , 17 ' 18 and magnetic fielcU 7 - are also in- 
fluencing Ng . The number of vortices inside a hole can 
still increase when the magnetic field further increases 
after the transition^ 3 - Therefore, the transition between 
region-I and region-II is mainly determined by the satu- 
ration number Ng, which is temperature, geometry and 
magnetic field dependent. 

When the temperature is below 8.640 K in Fig.gfa), we 
notice a decrease of resistance maxima in a wide range of 
magnetic field. This abnormal effect is related to the sur- 
face superconductivity around hole edges and it reaches 
its maximum near H C 2 (corresponding to Hjj)J^ Using 
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Figure 4. (color online) Resistance as a function of reduced 
magnetic flux / for samples SI [upper panel (a)] and S2 [lower 
panel (b)], measured at several temperatures with fixed cur- 
rents I — 30 [iA and 200 /iA, respectively. From top to bot- 
tom, the corresponding temperatures of R(f) curves are 8.723 
K, 8.702 K, 8.680 K, 8.615 K and 8.540 K (partly shown in 
the figures). 



this feature, we have plotted Hu as a function of the 
reduced temperature, and obtained a linear temperature 
dependence, as shown in the inset of Fig. [SJb). It is 
found that Hr>(t) nearly coincides with the H C 2(t) for the 
Nb thin film. We use the thin film expression H* 2 (t) = 
<&o/[27r£(£) 2 ] to roughly estimate the upper critical field. 
The theoretical value H* 2 (T/T c0 = 0.986) = 361 Oe is 
close to H D = 352 Oe, where T = 8.640 K, and T c0 
= 8.763 K is obtained from a linear extrapolation of line 
Hn(t) to zero field. This suggests that Hd can be indeed 
identified as H C 2 of sample SI. The region between H c % 
and H C 2 is the region where the bulk sample is already in 
the normal state and only a superconducting sheath per- 
sists at the surface of the sample parallel to the applied 
field. At a relatively low temperature, as the edge states 
merge with wide strips containing a finite value of order 
parameter, the whole array reentrances to the supercon- 
ducting state. This has been found at T = 8.540 K for 
S2. Another possible reason for the abnormal behavior 
is lowering mobility of interstitial vortices by a 'caging 



effect ',22 

By comparing the R(H) curves of SI with S2, we can 
study the influence of geometry (aspect ratio) on the 
magnetoresistance oscillation. However, there is no much 
difference in the position of the fine structures in the 
low field regime or the position of dips in the high field 
regime, except the difference in the saturation number 
Ns and Hq. This implies that the oscillations in the 
R(H) curves are mainly influenced by the conductivity 
of the arrays and upper critical fields. According to the 
above study, we find that the transition boundaries Hg 
and Hd are both temperature and geometry dependent. 
It is worth to note that the concept of rectangular array 
should hold in a certain range of aspect ratio (or length 
of b). Too large aspect ratio turns the system to one di- 
mensional array of antidots, where the magnetoresistance 
behavior are dominated by the flux flow in the thin film. 
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C. Discussion 

The well-defined shape of dips enables us to accurately 
determine AH between two consecutive dips. FigureEJc) 
shows AH as a function of index number N for samples 

51 and S2, at T = 8.640 K and 8^5 K. The horizontal 
dashed lines (c and d) indicate AH for samples SI and 

52 in region-I. We find that AH = 20.8 ± 1 Oe for SI, 
and AH = 12.6 ± 0.9 Oe for C2. Two dashed lines (a 
and b) mark the theoretical values AH$l = $o/( 7r ^e//) 
= 37.9 ± 2 Oe for the singlc-loop-likc oscillations and 
AHsq = &o/a 2 =32.3 ± 1.7 Oe for square vortex lat- 
tice configuration^ The calculations are carried out with 
r e fj = 417 ± 10 nm and a lattice constant a = 800 ± 
20 nm. It is found that most data are larger than the 
level of AHsq and fall into the region of AHsl at higher 
field. The single-loop-like region marked with solid rect- 
angle in Fig. [5^c) are full of data points larger than 35 
Oe. The non-constant intervals are also consistent with 
the previous discussion. To further confirm our observa- 
tions, samples with smaller unit cells are made to obtain 
a larger AH in region-Ill. The results for samples Dl (400 
x 960, r h = 170 nm) are illustrated in Fig.EKd). For Dl, 
it is found that Ns = 2, which is smaller than that for 
SI, due to a smaller ru- AH values in the range of 272 
Oe < H < 402 Oe are close to the period of square lattice 
H Sq = $o/a 2 = 129.4 ± 6 Oe. In the range of H > 405 
Oe, AH is 143.6 Oe, which is close to the period calcu- 
lated for the single-loop-like oscillations, with ro = 216.5 
nm ~ r e ff = 230 nm. Similar results are also obtained 
in other samples (e.g. 400 x 655, 400 x 800 and 400 x 
1000). Therefore, we have observed a series of crossover 
behaviors for the oscillations of magnetoresistance, like 
the samples SI and S2. 

Finally, we compare our results with previous works. 
Firstly, in contrast to the square arrays of antidots, an 
extra transition (interstitial vortex state) between collec- 
tive oscillations and single-loop-like oscillations is found 
above the saturation field Hs- At higher fields, it in- 



Figure 5. (color online) (a) Temperature dependence of the 
saturation numbers Ns for sample SI and S2, contrasted with 
the theoretical values Nst(t) — r/i./2£(i). (b) Phase boundary 
for arrays of antidots and a reference Nb thin film. H c s is 
calculated as 1.695_ff C 2, where H C 2 is the experimental data 
for Nb thin film. H* 3 is the experimental data obtained from 
90%R n criterium in Fig-HJa). (c) AH values as a function of 
index number N for sample SI (A) and S2 (O), are taken at T 
= 8.640 K and 8.615 K, indicated by solid and open symbols, 
respectively. The horizontal dashed lines a to d mark the field 
values of 37.9, 32.3, 20.8 and 12.6 Oe. (d) R(H) curve for the 
sample Dl with rectangular array of antidots (400 nm x 960 
nm) measured at t — 0.983 with / = 50 /iA. 



dicates that the crossover from interstitial vortex state 
to the single-loop-like state is triggered by H c2 of the 
sample, offering another way to induce crossover behav- 
iors. Secondly, changes in the periodicity and the shape 
of the dips in magnetoresistance have been found in 
superconductors with magnetic dots, nonmagnetic dots 
and antidots i 21 ' 24 The dominant mechanisms discussed in 
those works are analyzed in terms of two possible models: 
the reconfiguration model and the multivortex model. 
Following these pictures, our results require that vortices 
form a denser state in high field region. Whereas, there 
are difficulties with the explanation of the phenomena, 
such as the periodic appearances of the fine structures in 
region-I, the non-hysteresis effects in region-Ill and the 
much large field intervals (42 Oe > AHs q , for SI) at 
fields higher than the upper critical field. On the other 
hand, from the viewpoint of dynamics of the vortex lat- 
tice ordering^ a monotonic increase in the intervals with 
increasing magnetic field is supposed to occur. However, 
this has not been found in our work. Thus, it seems to 
be inappropriate to use the London limit at such high 
temperatures and magnetic fields. We can discuss the 
results in the framework of the Ginzburg-Landau the- 
ory by considering the order parameter modulation. To 
clarify these effects, further direct imaging experiments, 
transport measurements and theoretical simulations in 



this type of arrays of antidots are necessary. 



IV. CONCLUSION 

In conclusion, we have investigated magnetoresistance 
of superconducting Nb thin films containing rectangular 
arrays of large antidots. The R{H) curve is divided into 
three regions by comparing the results with hysteresis 
measurements. At low magnetic fields, the system be- 
haves like a weak link wire network, giving rise to dips in 
R{H) at integral and fractional reduced magnetic flux be- 
low the saturation filed Hs- At the intermediate fields, 
the interstitial vortices form sublattice inducing larger 
magnetic field intervals in R(H) and a hysteresis effect in- 
dicates the existence of the barrier for the vortex motion. 
As soon as the magnetic field exceeds H c2 , superconduc- 
tivity nucleates near the edge of antidot in a layer with 
a width of £(T), resulting in a single- loop-like supercon- 
ducting state. In this state, the non-periodic oscillations 



and a fully reversible behavior arc found in the R{H) 
curves. The field intervals are determined from the effec- 
tive diameter of hole, which is comparable with the short 
side constant of the unit cell. The crossover fields (Hs 
and Hjj) between the three regions are both temperature 
and geometry dependent. 
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